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Abstract 

A lattice QCD calculation of the B meson decay constant is presented. In 
order to investigate the scaling violation associated with the heavy quark, 
parallel simulations are carried out employing both Wilson and the 0(a)- 
improved clover actions for the heavy quark. The discretization errors due 
to the large b quark mass are estimated in a systematic way with the aid 
of the non-relativistic interpretation approach of El-Khadra, Kronfeld and 
Mackenzie. As our best value from the quenched simulations at (3=5.9, 6.1 
and 6.3 we obtain /5=163ibl6 MeV and /bs=175±18 MeV in the continuum 
limit where the error includes both statistical and systematic uncertainties. 
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The B meson decay constant fs is a fundamental quantity needed to extract the Cabibbo- 
Kobayashi-Maskawa matrix element Vtd from experiments on — B^ mixing. For this 
reason lattice QCD calculations of fs have been pursued over several years, employing 
either relativistic or non-relativistic @ (including the static [Q) formulation for the b 
quark. 

While there are a number of advantages with the relativistic formulation , its basic 
problem for calculations of fs has lain in the difficulty of controlling systematic errors 
associated with heavy quark mass whose magnitude in lattice units exceeds unity for the 
b quark for a typical lattice spacing ^2 — 3 GeV accessible in current simulations. 
The formalism proposed in Ref. however, has shed a new light on this problem: these 
authors have pointed out that a Wilson-type lattice quark action for heavy quark can be 
reinterpreted as a non-relativistic Hamiltonian for an effective heavy quark field Q of the 
form, 



H = Q 
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where the effective heavy quark mass parameters rrii {i = 1,2,5, • ■ ■) are functions of the 
bare quark mass mq and the coupling constant. In contrast to the continuum where these 
parameters are equal (mi = m2 = ms = ■ ■ ■), they mutually differ by 0{amQ) at finite 
lattice spacing, which represents 0{amQ) errors of the original action in the framework of 
the effective Hamiltonian (|1|). These parameters, however, are calculable in perturbation 
theory, and effects of 0{amQ) errors on fs can be systematically analyzed. In particular, 
we observe that errors of 0{{m2/mB — l)AQCD/niQ) for the Wilson action (m^ 7^ 1712) is 
reduced to O^asA^Qco/niQ, K^f^j^/mq) for the 0(a)-improved clover action p|, for which 
ruB = m2 holds at the tree level. 

In this article we report on a calculation of the B meson decay constant in quenched 
lattice QCD with the relativistic formalism employing this "non-relativistic interpretation" . 
In order to study ©(amg) systematic errors, we carry out a parallel set of simulations using 
both Wilson and clover quark actions over a wide range of heavy quark mass and lattice 
spacing. 

The parameters of our simulations are listed in Table |. The standard plaquette action is 
employed for gluons. For the clover coefficient we use the tadpole-modified |^ one-loop value 
Csw = l/^o[l + 0.199ay(l/a)] where Uq = P^^^ with P the average plaquette. The lattice 
size is chosen so that the physical size is approximately kept at L ^2 fm. Seven values of 
the heavy quark hopping parameter are taken to cover the charm and bottom quark masses, 
and four values for light quark in a range 0.4ms — 1.4ms with ms the strange quark mass. 
The simulations have been carried out on the Fujitsu VPP500/80 supercomputer at KEK. 

We extract the heavy-light decay constant fp from the correlators of the axial vector 
current and the pseudoscalar density P given by {Ai{t)P{Q)) and {P{t)P{0)) . In order 
to reduce statistical errors of the correlators, which rapidly increase as amq increases, we 
employ the smeared pseudoscalar density P^{x) = J2r4>{\^)Q{^ + ''")l5l{x) on the gluon 
configurations fixed to the Coulomb gauge. The smearing function 0(|r|) is obtained by 
measuring the wave function of the pseudoscalar meson for each set of heavy and light 
quark masses. As a result we are able to isolate the ground state signal from a small time 
separation of t ~ 0.8 fm. 
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We adopt for the heavy-hght axial vector current q^fi^y^Q the one-loop renormalization 
factor Z^(amQ) newly calculated with full inclusion of the heavy quark mass dependence 
I0| . The calculation is available for both Wilson and clover actions, and it confirms the 



results of Refs. pT| , p!2| made earlier for the Wilson action. We find that effects of finite amq 
are non- negligible: with Z^(amQ) evaluated with the coupling constant Q;y(l/a), Jb for the 
Wilson action is reduced by 5% (at (3 = 5.9) to 2% (at j3 = 6.3) compared to the value 
obtained with the Za factor with the mass dependence ignored, as employed in the previous 
studies. For the clover action the finite amq effect works in the opposite direction with a 
similar magnitude. 

We remark that the field Q is related to the original field \E' through 

Q = e-'^^/^[l + d,j-D]^, (2) 

where di is a known function of amq [Q. The KLM factor [^] e"™^/^ is evaluated including 
the mg-dependent one-loop correction |jlO|]. We ignore the di^ ■ D term, since its corrections 



to fs is expected to be at most 1-2% due to a small value of c/i 0.1. 

A non-trivial issue in lattice studies of heavy-light mesons is how to define their masses, 
since the pole mass directly measurable from meson propagators with zero spatial momentum 
suffers from large 0{amq) errors. A possible choice is the kinetic mass defined by an 
expansion of the energy-momentum dispersion relation of the meson, 

^meson (p) = "^pole + 7) ^ ^ (p *^ ) . (3) 

The kinetic mass mkin, however, receives corrections from 0(p ^) terms in ([1|) which are 
uncontrolled and hence suffer from a large 0{amq) effect [|13|. This leads to a pathology 
that the h quark mass cannot be determined consistently from heavy-light and heavy-heavy 
mesons 0,0. 

An alternative choice may be to define a "kinetic mass" by taking the pole mass for a 
meson corrected by the difference of the kinetic and pole masses of the heavy quark m2 — mi 



rrikin = rripoie + ("^2 - rrii) (4) 

This choice is motivated by the expectation that the binding energy of a heavy-light meson 
becomes independent of the heavy quark mass in the non-relativistic limit and (m2 — mi) 
should thus represent the difference between kinetic and pole masses of the meson. We find 
that the meson mass calculated in this way does not suffer from the pathology observed with 
mkin defined by (pi). We adopt this definition in our analyses using the one-loop perturbative 



result [|10| for m2 — mi. 

Let us now present our results. We plot $(mp) = {as{mp) / asims))'^^^'' fpy/mp in Fig. |l| 
as a function of the inverse heavy-light meson mass mp for both Wilson (open symbols) and 
clover (filled symbols) actions. The light quark mass is linearly extrapolated to the chiral 
limit, and as(/i) is calculated with the standard 2-loop definition where we employ the value 
AqcD = 295 MeV estimated from the ay coupling using the plaquette average 

There is an ambiguity in practice as to what mass scale is to be adopted to represent a 
quantity that has mass dimension. We prefer to use the scale that does not depend on the 
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quark sector to facilitate a direct comparison of the 0{amQ) errors with the two different 
quark actions on the common gauge configurations. Hence our natural choice is the string 
tension a, and the ordinate is normalized by a^/^ and the abscissa by in Fig. |1|, where 
we employ the string tension of Ref. [|l^]. Vertical lines indicate the positions of the B and 
D mesons if one uses a phenomenological value y/a = 427 MeV [|17|. Data points plotted at 
l/mp = are the static results ||l8l, to which our data seem to converge towards the heavy 
quark mass limit. We observe that the Wilson results exhibit a small increase as the lattice 
spacing decreases, while the clover points at three values of f3 lie almost on a single curve. 

An improved scaling behavior with the clover action is more clearly seen in Fig. ^ where 
we present the continuum extrapolation of J'b and fo- Compared to scaling violation of 
11-5% in our range of lattice spacing a^^ ^ 1.6-3 GeV for the Wilson case, the clover 
data show a significantly reduced variation of 4-2% over the same range of lattice spacing. 
These magnitudes are common to Jb and fn- The continuum values obtained by a linear 
extrapolation agree within the statistical error of about 5% between the two actions. 

This agreement, however, does not necessarily mean that systematic 0{amQ) errors are 
all removed by the continuum extrapolation. Let us discuss this point for the Wilson action 
for the heavy quark. According to the non-relativistic Hamiltonian (|l]), the size of the leading 
O^amo) error in is 0{{cb — 1)Aqcd/^q) where cb = m2/mB- The tree level value cb = 
1/(1 + sinhmia) is plotted in Fig. ^ as a function of = e"^^°' sinhmia/(l + sinh rriia). 
For ~ 2.9 — 1.5, corresponding to the b quark at /5 = 5.9 — 6.3, \cb — 1| ~ 0.7 — 0.5, 
and hence we expect an error of 0(4 — 3%) in J'b at our simulation points. If we linearly 
extrapolate cb to the continuum limit = 0, \cb — 1| decreases to 0.4, which implies an 
0(3%) error left unremoved. For the D meson, the value of \cb — 1| is smaller {\cb — 1| ~ 
0.4 — 0.3 for the charm quark at m2a ~ 0.9 — 0.5) and decreases faster, extrapolating to 
|cb — 1| ~ 0.2 at m2a = 0. Thus, 0{amQ) errors of 0(7 — 5%) for fc, at our simulation 
points reduces to 0(3%) in the continuum limit. This consideration indicates that the use of 
non-relativistic Hamiltonian inherently leaves a mg-dependent systematic error that cannot 
be removed by a linear extrapolation. We estimate that it is of the order of 3% for Jb and 
/d in the continuum. 

We need to consider two more sources of systematic errors, which can in principle be 
removed by the extrapolation procedure if the simulation is made at high precision but in 
practice are not removed from our results due to the insufficient statistics. One of them 
is mg-independent scaling violation, which is O^aAgcD) for the Wilson action. We take 
the value aAqcn ~ 10% at our smallest lattice spacing a^^ ^3 GeV as an estimate of 
O^aAQcn) scaling violation effects. The other is the 0{ay) uncertainty due to the use of 
one-loop value for Za, which is 0(4%) with ay(l/a) ~ 0.2 at ^ 3 GeV. Therefore, we 
expect a systematic error of the order of 10% in our results for the decay constant obtained 
with the continuum extrapolation. 

This error analysis gives us some insight about the origin of the scaling violation observed 
in Fig. 1^. We can conclude that the dominant part of the lattice spacing dependence comes 
from the mg-independent oAqcd effect, since the mg-dependent errors 0{{cb — 1)Aqcd/^q) 
diminishes only little towards the continuum limit and hence contributes little to the slope as 
a function of a. This leads us to expect that Jb and Jd exhibit a similar slope as a function 
of a, as we indeed observe in the figure. The size of scaling violation actually observed is 
within a factor of two from our estimate. 
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For the clover action the mg-dependent errors are reduced to O^asAQCD/fng) and 
0{{Aqcd/^q)'^)- We estimate them to be 0(1%). The scahng violation error is O^asaAqcD) 
and 0{a'^AQ(-<^) which are of the order of 2%. Taking account of the 0{ay) error from Za 
and that arising from the field rotation term diKqcD/^Q ~ 0(2%) in (|^), which is ignored in 
the present calculation, we expect systematic errors of order 5% for the decay constant from 
the clover action. The mg-independent scaling violation also dominates the a dependence 
of the decay constant, the contribution of mg-dependent errors being very small. 

We now examine the question of how to set the physical scale of lattice spacing to 
calculate the decay constant. The most common in the literature is to use either p meson 
mass TUp or pion decay constant /,r to determine the lattice scale. In Fig. ^ we give the 
ratio of the lattice scale obtained with these quantities to that with the string tension. For 
the clover action the 0(a)-imp roved axial vector current A/^ + CAddiP is used to measure 
with the one-loop value for the coefficient ca 0. 

As expected, the slope of the ratio is much more gentle for the clover action compared to 
that for the Wilson action. The values in the continuum limit obtained by a linear extrapo- 
lation show a significant scatter, and the continuum limits of the ratio with the two actions 
disagree at the level of 5-10%. We ascribe this discrepancy mainly to smaller statistics of 
our Wilson simulation, and a resulting uncertainty in the continuum extrapolation. 

The continuum value of the ratio need not be equal to unity in the quenched approxi- 
mation; the disagreement may represent the systematic error due to quenching. Separating 
the quenching error from statistical and extrapolation uncertainties, however, is not possible 
with our present statistical accuracy. We then take the dispersion of the ratio in Fig. H as 
an uncertainty of the scale including the quenching error. We estimate it to be 10% for the 
Wilson action and 5% for the clover case. 

We present our results for the physical value of the decay constant in Table Here we 
set the scale using the p meson mass. To obtain the ratio fp/rrip in the continuum limit, we 
combine the continuum values of fpj \fo and ^fojmp obtained by a linear extrapolation as 
given in Figs, ^and ^ A direct continuum extrapolation of fp/mp yields consistent results. 
The errors quoted in the parentheses are, in the order given, statistical, systematic and scale 
errors. 

We take the result from the clover action to be our best estimate primarily because 
the uncertainty from scaling violation is smaller, but also because our statistical ensemble 
is larger for this case. Combining errors by quadrature we obtain /^=163±16 MeV and 
/b^=175±18 MeV for the B meson decay constants. For the D meson we obtain fn = 184±17 
MeV and /d, = 203 ± 19 MeV. 

We have shown in this article that B meson decay constant within a 10% accuracy can 
be obtained with the 0(a)-improved clover quark action in current lattice simulations at 
~ 1.6 — 3 GeV. The systematic error associated with the heavy quark is no longer the 
dominant source of uncertainty. The uncertainties in the lattice scale turns out to be more 
important in the present simulation. Time-consuming full QCD simulations are perhaps 
indispensable to go beyond the presently achieved accuracy in view of the fact that the scale 
uncertainty involving the quenching error will be the largest source of the uncertainty in the 
calculation of the heavy-light decay constant. 

This work is supported by the Supercomputer Project (No. 97-15) of High Energy Accel- 
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FIG. 1. as a function of 1/mp normalized by string tension a. Filled symbols represent 
results with the clover action and open symbols with the Wilson action. Circles, squares and 
diamonds correspond to results at /3=5.9, 6.1 and 6.3, respectively. Points at 1/mp = are static 
results [111] at the same set of /3. 
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FIG. 2. Continuum extrapolation of fp (circles) and fo (squares). Filled symbols represent 
results with the clover action and open symbols with the Wilson action. 
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FIG. 3. Tree-lcvcl evaluation of the coupling cb of the chromomagnetic interaction term in the 
non-relativistic effective Hamiltonian. Solid and dashed lines correspond to the clover and Wilson 
actions, respectively. 
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FIG. 4. Ratio of lattice scale obtained from rup (circles) and from /^r (squares) to that from 
the string tension. Filled symbols represent results with the clover action and open symbols with 
the Wilson action. 



9 



TABLES 



TABLE L Simulation parameters. The lattice scale quoted is estimated from mp=770 MeV. 



action 


13 


5.9 


6.1 


6.3 




size 


16^x40 


24=^x64 


32^x80 


Wilson 


Aconf 


150 


100 


100 




1/a (GeV) 


2.03(3) 


2.65(4) 


3.31(6) 


clover 


Aconf 


540 


200 


166 






1.580 


1.525 


1.484 




1/a (GeV) 


1.64(2) 


2.29(4) 


3.02(5) 



TABLE II. Results for the decay constant in MeV unit. 





AMlson 


clover 


fB 


140(11)(15)(24) 


163(9)(8)(11) 




159(10)(17)(27) 


175(9)(9)(13) 


fo 


163(13)(18)(28) 


184(9)(9)(12) 


Id, 


180(11)(20)(31) 


203(9) (10) (14) 
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